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ABSTRACT 

We present spectroscopic confirmation of a z = 0.99 galaxy cluster discovered using data from the 
Wide-field Infrared Survey Explorer (WISE). This is the first z ~ 1 cluster candidate from the Massive 
Distant Clusters of WISE Survey (MaDCoWS) to be confirmed. It was selected as an overdensity of 
probable z > 1 sources using a combination of WISE and SDSS-DR8 photometric catalogs. Deeper 
follow-up imaging data from Subaru and WIYN reveal the cluster to be a rich system of galaxies, 
and multi-object spectroscopic observations from Keck confirm five cluster members at z = 0.99. The 
detection and confirmation of this cluster represents a first step towards constructing a uniformly- 
selected sample of distant, high- mass galaxy clusters over the full extragalactic sky using WISE data. 

Subject headings: galaxies: clusters: individual (MOO J2342. 0+1301) — galaxies: distances and 
redshifts — galaxies: evolution 



1. INTRODUCTION 

Historically, clusters of galaxies have been used as pow- 
erful probes of cosmology and galaxy evolution, provid- 
ing such landmark result s as the first evidence for the 
existence of dark matter (Zwicky 1937), demonstration 
of the impor tance of environment m galaxy evolution 
( Dresslerjl980[), and direct proof of the existence of dar k 



matter ( |(Jlowe et al 
The unique leverage provi 



2004 [2006 
•roviaed by 



Bradac et al.||2006[ ). 
galaxy clusters comes 
primarily from their extreme mass (M > 10 14 M ) and 
late-time growth that continues up to the present epoch. 

Large-area surveys afford the opportunity to identify 
well-defined samples of the most massive, rarest galaxy 
clusters. The ROSAT All- Sky Survey, for example, has 
yielded several notable catalogs of massive X-ray se- 
lected galaxy clusters to moderate redshifts (e.g. BCS 
at z < 0.3 and MACS at z < 0.7; |Ebeling et aTlpMr 



2001), while the Sloan Digital Sky Survey (SDSS) has 



produced large catalogs of nearby clusters covering a 



wider range of cluster masses (e.g. |Koester et al. 2007 



0.1 < z < 0.3). The Planck mission also provides an 
all-sky catalog of very massive galaxy clusters at z < 1 
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(Planck Collaboration et al. 2011a|b ), while the South 
Pole Telescope (SPT) and Atacama Cosmology Tele- 
scope (ACT) provide complementary samples of high- 
mass clusters reaching to z > 1 for several thousand 



square degrees ([W illiams on^ et al.||2011| |Marriage et al 
2011 Reichardt et al.||2012| ). However, there currently 
exist no surveys capable of identifying massive clusters 
at z > 1 over the full extragalactic sky. 

The most massive clusters in this redshift regime are 
of particular interest given the recent vigorous debate 
about whether the few known massive clusters at z > 1 
are co nsistent with Gaussian primordial density fluctua- 



tions ([Cayon et al.|2Qlll|Hoyle et al.|20Tfl|Enqvist et al 
2011 Iwilliamson et al.p011| |Jee et al.||20ill |Hotchkis' 



2011 

smal 



nliamson et al.||2Ull| |Jee et al.||2011| |Hotchkiss 
but a definitive answer remains elusive due to 
number statistics. Moreover, the recent discovery 
of strong lensing by a galaxy cluster at z = 1.75 (|Stanford| 
et al.||2Q12| |Brodwin et~aLl|2Q12} |Gonzalez et al.l|2012p 
revives the question of whether the frequency of giant 
arcs behind clusters at high redshift is consistent with 
ACDM, a question which can only be comprehensively- 
addressed with a statistical sample of massive clusters in 
this regime. 

The various wavelength regimes and selection tech- 
niques used to construct galaxy cluster samples each offer 
unique advantages and disadvantages. Catalogs selected 
using the Sunyaev-Zel'dovich (SZ) effect provide high pu- 
rity samples that are nearly mass-limited due to the de- 
tection signal's weak redshift dependence. These features 
make SZ surveys ideal for cosmological tests based upon 
evolution of the cluster mass function. However, the cur- 
rent generation of SZ catalogs are limited to z < 1.3, 
and at z > 1 to M200 ^,3x 10 14 and a few thousand 



ana ai z ^ 1 10 1V1200 ° x iu ivip aim a lew inousana 
squar e degrees (e.g. |Marriage et al.|2011 Reichardt et al. 
2012). Similar to trie SZ, surveys at X-ray wavelengths 



provide an estimate of cluster mass immediately from de- 
tection. Current X-ray surveys are able to probe further 
down the mass function at high-redshift. Recent searches 
with XMM- Newton have proven sensitive to clusters out 
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Fig. 1. — WISE Wl— W2 vs. z for various simulated stellar pop- 
ulations. "Burst" models passively evolve after a single starburst 
lasting 0.1 Gyr. "Exp" models have a star formation history of 
the f orm e~*/ r with r = 10 G yr. Models labeled "BC03" come 
from |Bruzual fe Chariot] ( |2003| and have a[Chabrier (2003) initial 
mass function, while those labeled "M05 come rrom Maraston 
( 2005 ) and have a |Kroupa| ( |2001| ) initial mass function. The for- 
mation redshift is given by Zf. Trie dashed horizontal line denotes 
the W1-W2 cut value used to select MOO J2342. 0+1301. Mod- 
els we re generated using the EzGal package ( Manc one fc Gonzalez| 
|2012) . 



to z > 1.6 for systems with M200 > 10 14 M fi flFinoguenov 



et al. 120071 IPapovidi et al.|2QlQ|[Tanaka et al.|201U| ^'ass- 
bender et al.||201lj ). The key factor limiting the current 
li< 



X-ray surveys is the area surveyed, roughly an order of 
magnitude less than the SZ programs. Complementary 
to the SZ and X-ray surveys, cluster searches with Spitzer 
based upon detecting galaxy overdensities have, in re- 
cent years, provided the greatest reach in redshift and 
mass sensitivity at high-r edshift, extending to z > 2 and 
5 x 10 13 (e.g. Eise nhardt et al. 2008; G obat 
t al.||2U12[|Zeimann et al.||2012| ). 
Infrared Survey Explorer mission 



M 

let al.||2011| 
1 — rrn-L — rrr? 



The 



Stanfor d et 
Wide- field 



(WISE; IWright et al]j2010l ), which covers the full sky 
at 3.4, 4.6, 12, and 22 //m, offers the potential to find 
massive clusters at z > 1 over the full extragalactic sky. 
In this paper we present the first distant galaxy clus- 
ter discovered using WISE data, MOO J2342.0+1301 at 
z = 0.99. The discovery of MOO J2342.0+1301 repre- 
sents a first step towards constructing an all-sky sample 
a t 1 ^ z i$ 1-4 as part of the Massive Distant Clusters of 
WISE Survey (MaDCoWS). [^] 

2. WISE SELECTION OF CLUSTER CANDIDATES 

2.1. WISE Data 

The key aspect of the WISE mission enabling detec- 
tion of high redshift galaxy clusters is the excellent pho- 
tometric sensitivity at 3.4 and 4.6 fim (Wl and W2, re- 
spectively). The WISE All-Sky Data Release achieves 
5<j sensitivity limits better than 0.07 and 0.1 mJy in un- 
confused regions in Wl and W20 These sensitivities 

10 Sources in the MaDCoWS catalog are designated as MaD- 
CoWS Overdense Objects, which is the origin of the MOO target 
designat ion. 

11 [trttp : //wise g . ipac . caltech . edu/docs/release/al lsky/ 
expsup/ sec6_3a . html 



correspond to areas near the ecliptic plane and represent 
the typical minimum depth for the All-Sky survey. The 
effective exposure time of the observations vary signif- 
icantly with ecliptic latitude due to the survey design, 
increasing from twelve 7.7s Wl exposures at the ecliptic 
plane to more than a hundred near the ecliptic poles. The 
WISE All-Sky Source Catalog is based on this imaging 
data tiled into 18,240 Atlas Images of ~ 2.4 deg 2 , from 
which approximately 5.6 x 10 8 sources were extracted at a 
5a detection threshold. The WISE All-Sky Release data 
products, including both the Source Catalog and Image 
Atlas, became public on 2012 March 14 and can be ac- 
cessed via the NASA/IPAC Infrared Science Archive p 2 "] 

2.2. Sample Definition 

Full details of the MaDCoWS galaxy cluster search 
method will be presented in a future paper, but we pro- 
vide an overview of the essential features here. The 
search method that selected MOO J2342. 0+1301 was 
based on the W1-W2 color, following the approach of 
Papovich ( 2008[ ) , who discovere d a z = 1.62 galaxy clus - 



ter using Spitzer /IRAC data (Papovich et al. 2010). 
This method takes advantage oi the tact that galaxy 
colors in the observed- frame 3 — 5 /am regime become 
monotonically-redder between 0.75 < z < 1.75, and are 
largely insensitive to variations in star formation history 
(see Figure [T]). 

After cleaning the WISE catalog of flagged sources 
and compensating for spatial variations in depth, in 
the northern hemisphere we matched the WISE extrac 
tions to th e SDSS-DR8 photometric catalog (Aihara 
et al. 2011). For the preliminary search used to iden- 
tity MOO J 2342.0+1301, we then rejected sources using 
a combination of color and optical magnitude cuts, ex- 
cluding all objects for which W1-W2 < 0.3, z-Wl < 5, 
or i < 21. These cuts effectively remove the bulk of the 
foreground galaxy population at z < 1. The remaining 
sources, consisting predominantly of high-redshift galax- 
ies, were then binned into ~ 10° x 10° overlapping 2- 
dimensional density maps with a resolution of 15" pix -1 . 
These maps were smoothed with a Gaussian-difference 
wavelet kernel tuned to enhance structures on scales of 
~ 3 7 , corresponding to a physical size of ~ 1.4 Mpc at 
z ~ 1. From these maps we identified the most statis- 
tically significant overdensities for further investigation. 
MOO J2342.0+1301 was selected from - 10,000 deg 2 of 
WISE/SDSS-DR8 overlap and was ranked very highly 
by both raw overdensity and visual assessments by the 
MaDCoWS team. Figure [2] shows the cluster in the four 
WISE bands, centered on the location of the brightest 
pixel in the associated wavelet map peak. 

3. FOLLOW-UP OBSERVATIONS 

3.1. Subaru Imaging 

As part of the MaDCoWS follow-up program, 
MOO J2342.0+1301 was imaged u sing MOIRCS 
( |Ichikawa et a l. 2006; Suzuki et al.||2008| ) on Subaru in J 
and K s on UT 2011 September 18 and 19, respectively. 
A total of 20 minutes of integration time was obtained in 
both bands, with J observed in clear conditions and 0V6 
seeing and K s observed in slightly hazy conditions with 
0'/8 seeing. 



http : / / wise2 . ipac . caltech . edu/ docs/ release/ allsky 
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Fig. 2. — Top: 5' x 5' images in the four WISE bands, centered at 23h42m05.6s +13d01m29s. The cluster center is taken from the 
brightest pixel in its associated wavelet map peak. Bottom Left: 5' x 5' WISE-W1 image of MOO J2342. 0+1301. The white circles 
represent sources identified as probable high-redshift galaxies based upon color and magnitude cuts (see text). The circles have a diameter 
of 6.1 /x , comparable to the full width half at maximum of the WISE bands Wl and W2 point spread function. The box denotes a 3' x 3 ; 
region centered on the detection position. Bottom Right: Subaru-MOIRCS K s image of MOO J2342. 0+1301. The box and circles are 
identical to those in the left panel. Many of the sources in the WISE imaging are blends which resolve into multiple galaxies in the deeper 
follow-up imaging. 



The MOIRCS data were reduced using a set of IRAF 
script s developed by one of the authors ( jTanaka et al. 
2011| ). The reduction procedure follows the standard 
methodology for the reduction of IR imaging of faint 
sources. After flat fielding, sky subtraction, and cos- 
mic ray removal, the individual frames were corrected to 
remove geometric distortions, using a solution provided 
by the Subaru Observatory. These corrected frames were 
then registered and coadded for each detector and filter. 
Astrometric calibration and m osaicing of th e chips were 
performed using th e SCAMP ( |Bertin||20Q6| > and SWarp 
QBertin et al.|2002| software packages, using the Two Mi- 
cron All Sky Survey as the as trometric reference frame 
(2MASS; [Skrutskie et aLp QQ6). The data were also pho- 
tometrically calibrated using 2MASS. The K s image is 
shown in Figure [2j 



from the Mini- 
on WIYN dur- 
, A total of 60 



3.2. WIYN Imaging 

We obtained Sloa n i-band imagini 
Mosaic instrument (Saha et al. 2000| 
ing the night of UT 2011 INovember 
minutes of exposure was obtained in 0'/9 seeing using 
three dithered 20 minute integrations. The data were 
processed using standard reduction procedures in IRAF 
and astrometrically calibrated to the USNO-B1.0 catalog 
QMonet et aT]|2003| > using SCAMP and SWarp. Photo- 
metric calibration of final, stacked images was performed 



using SDSS-DR8 stars in the field. A color composite of 
the z, J and K s images is shown in Figure [3j 

3.3. Keck Spectroscopy 

Based on the Subaru imaging results, a slit-mask was 
designed to be used with the Low- R esolution Imaging 
Spectrometer (LRIS; |Oke et al |[l995| ) at the Keck Ob- 
servatory to obtain redshilts of potential cluster galax- 
ies. Sources for the mask were selected using J — K s 
colors, K s magnitudes, and distance from the nomi- 
nal cluster center, with final selection also accounting 
for slit-positioning constraints. Spectra were obtained 
on UT 2011 October 22 using l'/l x 10" slitlets, the 
G400/8500 grating on the red side, the D680 dichroic, 
and G300/5000 grism on the blue side. Four 1200 s ex- 
posures were obtained in mostly clear conditions with 
0'/6 seeing. 

The LRIS spectra were split into slitlets which were 
separately reduced using standard long-slit procedures 
in IRAF. The relative spectral response was calibrated 
via longslit observations of Wolf 1346 and Hiltner 600. 
Redshifts were determined by visual inspection using 
prominent spectral features including the 4000 Angstrom 
break, the Ca H+K absorption lines, and the [O II] A3727 
emission line. 

4. RESULTS 
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Fig. 3. Left: A WIYN-z and Subaru JK S 3-color image of the central 3' x 3' region of MOO J2342. 0+1301. Boxes indicate 
spectroscopically confirmed cluster members. The diamond denotes an object with a redshift that is not at the cluster redshift. Of the 
six galaxies targeted in the LRIS mask which are within 500 kpc of the nominal cluster center, five were confirmed to be cluster members. 
Right: LRIS spectra of the five member galaxies. The rest-frame wavelength at z = 0.99 is shown along the top. The vertical dotted lines 
mark the positions of the following features (from left to right): [O Il]A3727, Ca H+K, H5, and the G-band. 



4.1. Redshifts 

The Keck spectroscopy yielded high-quality redshifts 
for 12 galaxies. Geometrical constraints on the slit place- 
ment allowed observations of only six galaxies within a 
projected distance of 500 kpc from the nominal cluster 
center. Five were confirmed to be cluster members, with 
redshifts lying within Az = 0.010, or ±750 km s _1 , of 
z = 0.987, the mean cluster redshift. In Table [I] we 
provide coordinates, redshifts, identifying features, and 
a quality assessment for each redshift. All but one of 
the confirmed cluster members were identified based only 
on absorption lines and continuum breaks, rather than 
emission lines, a feature characteristic of the members of 
mature galaxy clusters. The spectra of the five confirmed 
cluster members are shown in Figure [3] 

4.2. Color-Magnitude Diagram 

We use imaging data from Subaru and WIYN to con- 
struct the J — K s color-magnitude diagram (CMD) for 
MOO J2342.0+1301. Ph otometry was obtain ed using 
SExtractor version 2.5.0 ( |Bertin fc Arnouts] |1996[ ) in 
dual-image mode with K s used lor source detection. We 
present the resulting CMD in Figure [4j Light grey dia- 
monds show all galaxies within a projected 500 kpc ra- 
dius of the cluster center that have i — J > 1.38, the 



color given by a Bruzual & Chariot (2003) exponential 
model with r = 10 Gyr, formed at Zf = 5 and observed at 
z = 0.99. Spectroscopically confirmed members and non- 
members are marked with black boxes (red in the online 
version) and black crosses, respectively. Additionally, we 
denote with large circles (blue in the online version) all 
i^ s -selected objects that lie within the footprint of the 
-6.1" WISE PSF for each WISE source that satisfies 
the color and magnitude cuts described in Section |2.2| 
(see also Figure [21. It is these WISE sources, the ma- 
jority of which are composed of multiple blended com- 
ponents, that contributed to the detection signal for the 
cluster. We include lines showing the expected J — K s 



color (horizontal) and apparent K s magnitude (vertical) 
of an L* galaxy at z = 0.99, using a |Bruzual fc Chariot 
( 2003[ ) sim ple stellar population model with solar metal- 
hcity and a |Chabrier| ( |2003| ) initial mass function, formed 
at Zf = 5! The red sequence in this cluster is well- 
populated, with both the spectroscopic members and 
galaxies corresponding to the WISE sources contribut- 
ing to detection of the cluster lying near the expected 
z 0.99 color. 

5. DISCUSSION 

The discovery of MOO J2342. 0+1301 represents the 
first result from MaDCoWS, whose goal is to conduct a 
wide-area search for such massive clusters at 1 < z < 1.4. 
This program is designed to be complementary to pre- 
vious, shallower all-sky cluster surveys, enabling future 
investigations of the most extreme clusters at this epoch. 
Forthcoming papers will present details for the full MaD- 
CoWS program, including the ongoing followup program 
and characterization of the cluster selection function. 

Finally, we note that this cluster was found using the 
WISE All-Sky Data Release, which includes only data 
obtained during the cryogenic mission. Data continued 
to be collected at 3.4 fim and 4.6 fim after exhaus- 
tion of the WISE cryogen supply as part of the NEO- 
WISE flMainzer et afl|2011|) survey of the asteroid belt. 
However, the asteroid search required only that individ- 
ual exposures be processed, so the multiple observations 
of inertially fixed sources from this phase of the mis- 
sion have not been combined. Fully processing and com- 
bining this post-cryogenic data with the cryogenic data 
in these bands, which has now been funded by NASA, 
would double the all-sky coverage from WISE in Wl and 
W2, improving the fidelity and redshift reach of WISE 
for discovering the most massive, distant galaxy clusters. 



The authors thank the anonymous referee whose com- 
ments improved the quality of the manuscript. This 
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Fig. 4. — (J — K s ) vs. K s color-magnitude diagram for MOO J2342. 0+1301. Gray diamonds denote galaxies that lie within a 500 kpc 
projected radius of the cluster center with i — J > 1.38. Black squares (red in the online version) are spectroscopically confirmed cluster 
members. Black "X" points represent non-member galaxies with spectroscopic redshifts. Large circles (blue in the online version) indicate 
all objects found within 3.05 /x of the WISE sources that contributed to the detection signal (see Figure The vertical and horizontal 
dashed lines show the expected K s magnitude and J — K s color for an L* galaxy formed at Zf = 5 and observed at z = 0.99. All magnitudes 
are on the Vega system. 



TABLE 1 

Spectroscopic Results 



a (J2000) 


5 (J2000) 


J 


K s 


z 


Quality 




Features 








Members 








23:42:04.80 


+13:00:50.0 


18.40 


16.75 


0.989 


A 


[O II], D4000, Ca H+K 


23:42:03.43 


+13:01:32.7 


18.66 


16.81 


0.993 


B 


D4000, Ca H+K 


23:42:05.67 


+13:01:24.2 


18.37 


16.69 


0.983 


A 


D4000, Ca H+K 


23:42:06.57 


+13:01:34.9 


19.33 


17.47 


0.987 


B 


D4000, Ca H+K 


23:42:07.56 


+13:01:56.8 


19.81 


17.97 


0.985 


A 


D4000, R8 








Non-Members 








23:42:03.93 


+13:00:53.2 


17.65 


16.00 


0.784 


B 


Ca H+K 


23:41:57.98 


+13:00:51.0 


19.17 


17.56 


1.043 


B 


[O II 




23:41:58.15 


+12:59:17.1 


18.66 


17.06 


0.777 


A 


Ca H+K, D4000 


23:42:08.86 


+13:02:47.9 


20.38 


18.47 


1.210 


B 


D4000 


23:41:59.87 


+12:59:14.9 


20.47 


18.94 


1.172 


A 


[O II 




23:41:59.28 


+12:58:09.6 


20.34 


18.48 


0.622 


A 


[O II 


, H/3 


23:42:00.65 


+12:58:16.2 


18.33 


16.87 


0.249 


A 


Ha, 


S II] 



Note. — Coordinates, redshifts, identifying features, and quality assessments of each 
redshift obtained from LRIS observations on Keck I. The typical redshift uncertainty is 
5z ~ 0.001. Quality flag "A" signifies a robust redshift determination, typically relying 
upon multiple emission or absorption features. Quality flag "B" signifies a redshift 
determination that is less certain, but still unambiguous with at least two features. 
Magnitudes are on the Vega system. 
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